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Abstract 

A theory of collisionless fluids is developed in a unified picture, 
where nonrotating (fii = 0,2 = ^3 = 0) figures with anisotropic 
(011 = 022 = C33) random velocity component distributions and ro- 
tating 7^ 7^ ^3) figures with isotropic [p\\ 7^ °22 7^ "733) ran- 
dom velocity component distributions, make adjoints configurations 
to the same system. R fluids are defined as ideal, self-gravitating 
fluids satisfying the virial theorem assumptions, in presence of sys- 
tematic rotation around each principal axis of inertia. To this aim, 
mean and rms angular velocities and mean and rms tangential veloc- 
ity components are expressed, by weighting on the moment of inertia 
and the mass, respectively. The figure rotation is defined as the mean 
angular velocity, weighted on the moment of inertia, with respect to 
a selected axis. The generalized tensor virial equations (Caimmi & 
Marmo 2005) are formulated for R fluids and further attention is de- 
voted to axisymmetric configurations where, for selected coordinate 
axes, a variation in figure rotation has to be counterbalanced by a 
variation in anisotropy excess and vice versa. A microscopical analy- 
sis of systematic and random motions is performed under a few general 
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hypotheses, by reversing the sign of tangential or axial velocity com- 
ponents of an assigned fraction of particles, leaving the distribution 
function and other parameters unchanged (Meza 2002). The appli- 
cation of the reversion process to tangential velocity components, is 
found to imply the conversion of random motion rotation kinetic en- 
ergy into systematic motion rotation kinetic energy. The application 
of the reversion process to axial velocity components, is found to im- 
ply the conversion of random motion translation kinetic energy into 
systematic motion translation kinetic energy, and the loss related to 
a change of reference frame is expressed in terms of systematic (imag- 
inary) motion rotation kinetic energy. A number of special situations 
are investigated with further detail. It is found that a R fluid always 
admits an adjoint configuration where figure rotation occurs around 
only one principal axis of inertia (R3 fluid), which implies that all 
the results related to R3 fluids (Caimmi 2007) may be extended to 
R fluids. Finally, a procedure is sketched for deriving the spin pa- 
rameter distribution (including imaginary rotation) from a sample of 
observed or simulated large-scale collisionless fluids i.e. galaxies and 
galaxy clusters. 

keywords - galaxies: clusters - galaxies: haloes - stars: stellar systems. 

1 Introduction 

Due to particle shocks, collisional fluids (e.g., stars, gas clouds) exhibit an 
isotropic stress tensor (cr^ = o\ 2 = o"| 3 ), where o^ v is the rms random ve- 
locity component on the axis, x p . The absence of particle shocks (leaving 
aside extreme situations, such as high-density galactic nuclei) makes a dif- 
ferent situation in collisionless fluids (e.g., galaxies, galaxy clusters), where 
the stress tensor is - in general - anisotropic {a\ Y ^ a\ 2 ^ 033)- The shape 
of the body is determined by systematic rotation, which is quantified by a 
spin parameter (null for nonrotating configurations), and/or by the difference 
between stress tensor diagonal components, or any equivalent anisotropy in- 
dicator, related to the rotation and an equatorial principal axis of inertia, 
respectively (null for configurations where the random velocity component 
distribution is isotropic). A description of collisionless fluids based on the 
equivalence of systematic and random motions with respect to the shape, 
appears to be highly rewarding and it would provide further insight on the 
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properties of stellar and galaxy systems. 

In an earlier attempt (Caimmi 1996) the stress tensor has been expressed 
as the sum of two terms, one related to a random (isotropic) velocity com- 
ponent distribution, and one other to anisotropic internal motions within 
the system. Further investigation has been devoted to the simplest situation 
where the system is made of two equal components, which are rotating at 
the same rate but in opposite sense. Then it has been recognized that the 
anisotropy excess may be related to real rotation, if the shape is flattened, 
and to imaginary rotation, if the shape is elongated, with respect to the 
rotation axis. 

A latter approach (Caimmi & Marmo 2005) has been restricted to homeoidally 
striated density profiles, for which the tensor virial equations were formulated 
and generalized to unrelaxed configurations. The kinetic-energy tensor has 
been expressed as the sum of two terms, one related to systematic rotation 
obeying an assigned law, and one other to the remaining motions e.g., ran- 
dom motions, streaming motions, radial motions. Finally, an expression of 
the spin parameter in terms of the anisotropy excess, has shown the role of 
systematic and remaining motions in flattening or elongating the shape. 

The above mentioned results have been improved and extended in sub- 
sequent work (Caimmi 2006, hereafter quoted as C06), where imaginary ro- 
tation has been related to negative anisotropy excess. Then sequences of 
configurations for which the generalized tensor virial equations hold, have 
been determined for homeoidally striated Jacobi ellipsoids including prolate 
shapes induced by imaginary rotation. The results from numerical simula- 
tions on the stability of rapidly rotating spherical configurations (Meza 2002) 
have been interpreted in the light of the theory. To this respect, the key argu- 
ment is that the reversion (from clockwise to counterclockwise or vice versa) 
of tangential velocity components related to an assigned fraction of partices, 
preserves the potential energy, the kinetic energy, and the distribution func- 
tion (Lynden-Bell 1960, 1962; Meza 2002). 

The study on homeoidally striated Jacobi ellipsoids has been extended to 
a more general class of bodies (R3 fluids) in a recent paper (Caimmi 2007, 
hereafter quoted as C07), where the contribution of radial and tangential 
velocity components on the equatorial plane was investigated with further 
detail. In addition, mean and rms (weighted on the moment of inertia) angu- 
lar velocity have been defined, and related to systematic and random motion 
tangential kinetic-energy tensor components, respectively. Also for R3 fluids, 
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it has been realized that the effect of (positive or negative) anisotropy excess 
is equivalent to additional (real or imaginary) figure rotation. 

The current attempt is aimed to extend the above mentioned results to a 
still more general class of bodies, R fluids, defined as ideal, self-gravitating, 
collisionless fluids where rotation occurs around each principal axis of iner- 
tia. It will be found that R fluids always admit an adjoint configuration 
where figure rotation occurs around a single principal axis, that is a R3 fluid. 
Accordingly, all the results which hold for R3 fluids may be extended to R 
fluids. 

The work is organized as follows. A number of basic definitions are pro- 
vided in Sect. 2, including the inertia tensor, the angular-velocity tensor, and 
the angular-momentum tensor. The generalized tensor virial equations for 
R fluids are formulated in Sect. 3. The microscopical analysis of systematic 
and random motions , for a collisionless fluid made of N identical parti- 
cles, is performed in Sect. 4, where a velocity component reversion process 
is defined, and a number of special situations are analysed in detail with 
respect to kinetic energy changes from random to systematic motions and 
vice versa. A procedure aimed to the derivation of the spin parameter distri- 
bution (including imaginary rotation) from an assigned sample of observed 
or simulated objects, is outlined in Sect. 5. Some concluding remarks are 
reported in Sect. 6. 



2 Angular- velocity and angular-momentum ten- 
sor 



In the special case of solid bodies, rotation is rigid and occurs around a single 
axis which, in turn, can remain fixed or change its direction. Accordingly, 
the angular momentum and the rotation kinetic energy read (e.g., Landau & 
Lifchits 1966, Chap. VI, §§31-33; hereafter quoted as LL66): 
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where J — (J\, J 2 , J3) is the angular-momentum vector, = (Qi, 0,2,^3) 
the angular-velocity vector, and I' the inertia tensor: 



Ks = / p(x 1 ,x 2 ,x 3 ) 
J s 



r=l 



d 3 S ; (3) 



related to the density profile, p, within the volume, S, being 5 rs the Kronecker 
symbol. The diagonal components of the inertia tensor, I' u , V 22 , I 33 , are the 
moments of inertia with respect to the axes, x±, x 2 , x 3 , respectively. 

In addition, the inertia tensor is symmetric of second rank, which implies 
the existence of a reference frame, (0' X 1 X 2 X 3 ), where the inertia tensor is 
diagonal (LL66, Chap. VI, §32): 

Z. = 8T.r„ ; ( 4 ) 

the coordinate axes coincide with the principal axes of inertia, and the di- 
agonal components define the principal moments of inertia. Accordingly, 
Eqs. (1) and (2) reduce to: 

J r = Y,I'„£l r ; r = 1,2,3 ; (5) 

s=l 

T rot = i (j^n? + I' 22 Vl\ + I> 33 Ql) ; (6) 

where, in addition (LL66, Chap. VI, §32): 

I' rr <I' ss + I' tt ; r = 1,2,3 ; s = 2,3,l ; t = 3, 1, 2 ; (7) 

with regard to the body under consideration. 

The inertia tensor has been defined in a different way, as (e.g., Chan- 
drasekhar 1969, Chap. 2, §9; Binney & Tremaine 1987, Chap. 4, §3): 

I„= f p( XuX2tX ,) XrX .d'S ; (8) 

J s 

and the combination of Eqs. (3) and (8) yields: 

3 

^rs firs ^ ' Irr 1-rs j (9) 
r=l 

I'rr = hs + ht ; T^S^t ; (10) 

r rs = -I rs ; r^s ; (11) 
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or: 



2I„ = r„ + r u -I'„ ; r^s^t ■ (12) 
Irs = -I'rs 5 r^s ; (13) 

which translates one formulation into the other (e.g., Bett et al. 2007). 

In the general case of (collisional or collisionless) fluids, rotation could 
occur different from solid-body, and around each principal axis of inertia. 
Let (0 X1X2X3) be a generic reference frame and (0'XiX 2 X3) a reference 
frame where the origin coincides with the centre of inertia, and the coordinate 
axes coincide with the principal axes of inertia. Let the coordinate axes, X ± , 
X 2 , X 3 , be defined as the principal axes. Let fli, Q 2 , Q3, be the angular- 
velocity vectors (to be specified later) related to the principal axes of inertia. 
Let £l rs be the component of the vector, £l r , on the coordinate The 
(3 x 3) tensor, Q rs , is defined as the angular- velocity tensor of the system 
under consideration, with respect to the reference frame, (0 x\ x 2 x 3 ). Let 
uJi, U2, <^3, be the angular- velocity vectors related to the coordinate axes, x±, 
x 2) x 3 . The following relation holds: 



u s = 

r=l 



J2^rs ; s=l,2,3 ; (14) 



and the angular- velocity tensor, u rs = Q rs , can formally be defined. Simi- 
larly, the (3 x 3) angular-momentum tensor is expressed as: 

J'rs = I'rs^rs ', (15) 

J' s = jZKs^rs ; (16) 
r=l 

where the inertia tensor is related to the reference frame, (0 x-y x 2 x 3 ). 

In the special case where the reference frame, (Oxi x 2 x 3 ), coincides with 
(0' XiX 2 X 3 ), then £l rs = 5 rs VL rr . Accordingly, Eqs. (14)-(16) reduce to: 

Us — &s — &ss — ^ss '■> (17) 
'^rs ^rs-^rr^rr i (1^) 

f s = j s = r ss n s ; (19) 

where I' tt = I rr + I ss , r ^ s ^ t, represents the moment of inertia with respect 
to the principal axis of inertia, Xt- From this point on, it shall be intended 
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that the origin coincides with the centre of inertia, and the coordinate axes 
coincide with the principal axes of inertia. 

The rotation kinetic-energy tensor is defined as: 



1 / 1/2 

(T TO i) rs — I rs Q r Q s —5 rs I rr Vl r , 



(20) 



where the diagonal components of the angular-velocity tensor are expressed 
as (C07): 



= Q r = — f 
IL Js 



Q r (x 1 ,X 2 ,X 3 ,t) 

w r = (x^ + x^ ■ 



n r (xi,x 2 ,x 3 ,t) WrP{xi,x 2 ,x 3 ,t)d 3 S; r ^ s ^ t; (21) 
V4 >r (x 1 ,x 2 ,x 3 ,t) 



W r 



(22) 
(23) 



and fl r (xi,x 2 , x 3 , t) is the mean value related to all the particles at the time, t, 
within the infinitesimal volume element, d 3 S = dxi dx 2 dx 3 , centred on the 
point, P(xi,x 2 ,x 3 ), v^ r is the tangential velocity component on the (Ox s x t ) 
principal plane, and the moment of inertia, I' rr = I ss + I u , reads: 



I'„= / w 2 r p(x 1 ,x 2) x 3) t)d?S ; 
^ s 



(24) 



as expected from the theorem of the mean, in connection with Eq. (21). 

Similarly, the mean square diagonal components of the angular-velocity 
tensor are expressed as: 



Vt r = (9? r ) = — \tt r (x 1 ,x 2 ,x 3 ,t)\ w r p(x 1 ,x 2 ,x 3 ,t)d S; r ^ s ^ t; (25) 

and the related variance reads: 

{^) 2 = W)-&f ; (26) 

as known from statistics. 

At this stage, it may be useful to extend and generalize the definition of 
figure rotation. 

Figure rotation. Given a R fluid, the figure rotation is defined as the mean 
angular velocity, weighted on the moment of inertia, with respect to a selected 
principal axis. 
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In terms of tangential velocity components, the counterparts of Eqs. (21), 
(25), and (26) read: 



= f ^t{xi,x 2 ,x 3 ,t) w r p(x 1 ,x 2 ,x 3 ,t)d 3 S ; r ^ s ^ t ; (27) 
(vl r ) = — J s [Q r (x 1 ,x 2 ,x 3 ,t)\ wjp(x 1 ,x 2 ,x 3 ,t)d 3 S; r ^ s ^ t; (28) 



(29) 



and the combination of Eqs. (25) and (28); (21) and (27); (26) and (29); 
yields: 



M(vl) = I'„(QJ) ; 
Mi^ r f = r rr .(a r r ; 



(30a) 
(30b) 
(30c) 



which relate tangential velocity components on the (0 x s x t ) principal plane, 
to angular velocity components on the x r principal axis. 



3 The generalized tensor virial equations for 
R fluids 

Let R fluids be defined as (collisional or collisionless) ideal self-gravitating 
fluids where figure rotation occurs around all the three principal axes of 
inertia. Let (0 £1X2X3) be a reference frame where the origin coincides with 
the centre of inertia, and the coordinate axes coincide with the principal axes 
of inertia. Then the mean radial velocity components must necessarily equal 
zero: 

*V = ; KJ = (a WrWr ) 2 ; (31) 

where v Wr is the radial velocity component on the (0 x s xt) principal plane, 
perpendicular to the x r principal axis. Let positive and negative radial ve- 
locity components be defined as directed outwards and inwards, respectively. 
The same holds for the mean tangential velocity components: 

Wr = ® ; W) = ^rf ; (32) 
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even in presence of systematic rotation. Let positive and negative tangential 
velocity components be defined as rotating counterclockwise and clockwise, 
respectively. 

The kinetic-energy tensor may be expressed as the sum of two contribu- 
tions: one, related to systematic motions, and one other, related to random 
motions (C07). The result is: 

Tk s k t — (T sys ) kskt + (T rdm ) kskt ; k = w,(j> ; (33) 

where the terms on the right-hand side, using Eqs. (30) and (31), can be 
expressed as: 

{T S ys)w a wt = ; (34a) 

1 2 

(T rdm ) WsWt = -5 st M (cr WsWs ) ; (34b) 



1 ~ 

(Tsys)^ = -5 st r ss (n s ) 2 ; (35a) 

keeping in mind that nondiagonal components are null in the case under 
discussion, only diagonal components shall be considered from this point on. 
The combination of Eqs. (26) and (35) yields: 

T^ r = \l' r M) ; (36) 

which depends on the density profile via the moment of inertia, I' rr , and 
the tangential velocity component distribution via the mean square angular 
velocity, (£1%), regardless from the fraction of systematic and random motions. 

In terms of the contributions related to the axial components of the 
kinetic-energy tensor, T ss and T u , Eqs. (33), (35), and (36) read: 

(T^M = l -I u (W) ; £ = s,t ; (37a) 

1 2 

[(T sys )4, r 4, r ]u = -Iu(£l r ) ; £ = s,t ; (37b) 
[(T rdm )^ r ]« = \hi [(fif) - (fi r ) 2 ] ; e = s,t ; (37c) 
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where Eq. (10) has been used. 

The invariance of a vector with respect to a change of the reference frame, 
implies the validity of the relations (C07): 



ivl ) + (vl 



ir, ' ■•'<■■ \ V s) + ( V t) j 



(c~MJ r «J r ) {P<f> r 4>r) 



+ {?*)' 



(38) 
(39) 
(40) 



where the velocity components on the x s and x t principal axes are labelled 
by the indices, s and t, respectively. 

The combination of Eqs. (26), (29), and (38)- (40) yields: 



{a WrWr f = (a ss ) 2 + {o tt f - |p [(n ?) _ (n 



which makes Eqs. (33) and (34) translate into: 

'(of) - (n r f 



l(T sys ) 



[(^rdm) 

•i = o ; 



\u = ^Ma 2 u 
e = s,t ; 



(41) 



£ = s,t; (42a) 
(42b) 



in terms of the contributions related to the axial components of the kinetic- 
energy tensor, T ss and T tt . 

The generalized tensor virial equations of the second order can be formu- 
lated, extending the procedure used for R3 fluids (C07). The result is: 



(Q s ) 2 + (Q t ) 2 + M( rr (T 2 + (E pot ) rr = ; 



°H + °22 + ^33 > 



(43) 
(44) 



C PP = ^^ = ^ ; P= 1,2,3 ; (45a) 

Cii + C 22 + C33 = ^ = ^ = C ; (45b) 

where mean angular velocity components on the x r principal axis are due 
to systematic rotation around x s and £1X6 S. (E pot ) rr is the self potential- 
energy tensor, ( rr may be conceived as generalized anisotropy parameters 
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(Caimmi & Marmo 2005; C06; C07), and T r d m is the effective random kinetic 
energy i.e. the right amount needed for an instantaneous configuration to 
satisfy the usual tensor virial equations of the second order, defined by the 
effective anisotropy parameters (C07): 

l PP =^f^ = C f ; p= 1,2,3 ; (46a) 

J- rdm *» 

Cn + C22 + C33 = 1 ; (46b) 

and the condition, £ = 1, or £ pp = ( pp , p = 1, 2, 3, makes Eqs. (43) reduce to 
their standard counterparts. To get further insight, a microscopical analysis 
is needed. 

In the special case of axisymmetric configurations, In = I22, (-Spot) 11 = 
(-Spot) 22, and the combination of the related tensor virial equations, expressed 
by Eq. (43), yields: 



I PP 



(Q q ) 2 -(Q P ) 2 = Ma 2 (( qq -( pp ) ; p=l,2 ; q = 2, 1 ; (47) 



where a figure rotation excess, [(£\) 2 — (^ P ) 2 ], is counterbalanced by an 
anisotropy excess, (( qq — ( pp ); in particular, a null figure rotation excess 
implies a null anisotropy excess and vice versa. Accordingly, a flattening 
on the (Ox p x r ) principal plane, induced by the figure rotation excess, has 
to be counterbalanced by an elongation on the x q principal axis, induced by 
the anisotropy excess, to yield an axisymmetric configuration with respect 
to the x r principal axis. 



4 Microscopical analysis of systematic and ran- 
dom motions 

Given a collisionless R fluid, let N be the total number of particles and 
m the mean particle mass in absence of mass segregation i.e. local and 
global mean particle mass coincide. For simplicity, the equivalent description 
(C07) involving iV identical particles of mass, to, shall be considered. With 
regard to the (0 x s x t ) principal plane, let v$ r be the tangential velocity 
component on the above mentioned plane. It is worth noting (Meza 2002) 
that the distribution function is independent of the sign of t^ r , and the 
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whole set of possible configurations are characterized by an equal amount of 
both kinetic and potential energy. Numerical simulations show that spherical 
systems, even if rapidly rotating, are dynamically stable after reversion of 
the tangential velocity component in an assigned fraction of particles (Meza 
2002). 

For sake of simplicity, let the initial configuration be nonrotating (v^ = 0) 
and with isotropic random velocity component distribution (£n = (22 = (33)- 
In the case under discussion of identical particles, — m, 1 < % < N, 
the centre of inertia velocity components, vc r , equal the related arithmetic 
means: 



N 



m 



V C r = 



1=1 



N 

i=i 



N 

E m 

i=i 



(0 



Nm 



(48) 



and the moments of inertia, I' rr , reduce to: 



r 



N 

i=i 



w 



m 



N 

E K i) 



1=1 



MR 2 Gr ; 



(49) 



where Rc r = [w r 2 ] 1 '' 2 is the curl radius with respect to the x r axis. 

The weighted mean, mean square, and rms tangential velocity compo- 
nents, expressed by Eqs. (27)-(29), read: 



N 



N 



Mh ^ ~ M 



E mw ^ 
i=i 

(<) = ^X> (0 



i=l 



E< ; 

i=l 

2 m 
M 



N 



— V \v {€) 

1 = 1 

fa&y = «) - (^) 2 = k^j 2 ^ 

and Eqs. (30) reduce to: 



= «) 



{Ptfrtrf = R Gr ( a fT r ^) I 



(50) 

(51) 
(52) 



(53a) 
(53b) 
(53c) 
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which relate weighted angular velocities around the to mean tangential 

velocity components on the (Ox s x t ) principal plane. 

At this stage, let the tangential velocity component of a fraction, n/N, 
of particles, be reversed in equal sense (from clockwise to counterclockwise 
or vice versa), according to the following assumptions. 

(i) Both the number, n, of particles where the tangential velocity component 

has been reversed, and the number, N — n, of particles which remain 
unchanged, are sufficiently large, 1 C n < iV, < n < Int(iV/2). 

(ii) The fraction, nk/Nk, of particles where the tangential velocity compo- 
nent has been reversed, within a generic volume element, Sk, is inde- 
pendent of the volume element, nk/Nk = n/N. 

(iii) The system is made of identical particles, — m, M — mN . 

(iv) After tangential velocity components have been reversed in rik particles 
within a generic volume element, Sk, on a total of Nk, a second set of 
7ifc particles (among the remaining Nk—rik) exists, where the tangential 
velocity component of any particle equals its counterpart belonging to 
the first set. 

In the following, the above process shall be quoted as "the reversion process" . 

Obviously, mean square tangential velocity components, (i>2), are left un- 
changed by the reversion process. On the contrary, mean tangential velocity 
components after the reversion process read: 



1 N i 
1 V- (i) _ 1 



2n N 



(54) 



1 i=2n+l 

where the first sum within brackets relates to particles where the reversion 
process has occurred and their counterparts with equal tangential velocity 
components, while the second sum comprises the remaining particles and 
necessarily equals the mean tangential velocity component before the occur- 
rence of the reversion process, which is null in the case under discussion. 
Accordingly, Eq. (54) reduces to: 

2n 

Wr = jf(v£)n ; (55a) 

1 2n -i n 

z ' i i=i i=i 
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keeping in mind that the first sum is performed on couples of particles with 
equal tangential velocity components. 

The validity of Eqs. (54) and (55) still maintains if tangential velocity 
components, v^ r , are replaced by axial velocity components, v r . The combi- 
nation of Eqs. (48) and (55a) yields: 

2n 

v Cr = — r = — (W)n ; (56) 

which is the velocity component of the centre of inertia with respect to the 
x r principal axis, after the reversion process. 

The total kinetic energy is left unchanged by the reversion process but, 
on the other hand, a fraction of random motion kinetic energy is turned into 
systematic motion kinetic energy. In the following, the reversion process shall 
be discussed with further details for a number of different situations. 



4.1 Tangential velocity component reversion 

Performing the reversion process on a given fraction of particles, n/N, with 
respect to tangential velocity components, implies the conversion of random 
(rotation) motion kinetic energy into systematic (rotation) motion kinetic 
energy, as: 

1 2n 

A(T rdm )^ r = -A(T sys )^ r = --M(v^f = --nm[(^) tt f ; (57) 

where the remaining parameters are left unchanged. 

The occurrence of the reversion process implies the following energy 
changes: 

2n 

T Td m -> T rdm - -nm[(^)„f ; (58) 

In 

(T rdm )^ r -> (T rdm ) 0r ^ - — nm[(v^) n } 2 ; (59) 

1 2n 

(T rdm )«^(T rdm )«---nm[(^) n ] 2 ; £ = s,t ; (60) 
In 

T sys -> + — nm[(v^) n f ; (61) 
In 

(TlysW -»• + J^ nm [(v^)n} 2 ; (62) 
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1 2n 

( T sy S )«^0 + -— nm(v^) 2 n ; £ = s,t ; (63) 

while the contributions from random radial motions along the equatorial 
plane, (T rdm ) WrWr , and the rotation axis, (T rdm ) rr , remain unchanged. 

With the system being relaxed, ( = 1, in the case under discussion, the 
generalized and effective anisotropy parameters, ( pp and ( pp , coincide with 
their counterparts related to the usual tensor virial equations, and Eqs. (45a) 
and (46a) take the explicit form (C06): 

t (l/3)T rdm -(2n/jV)(n/2)m[(ly;) n ] 2 

Qe = Tf, fey Mn w ^ 12 ' £ = s,t , (64a) 

T rdm - (2n/N)nm[{v (j>r ) n \ 



(1/3)71 



C = v ^"^ rdm ■ (64b) 

where the special case, n — 0, relates to the initial configuration, character- 
ized by isotropic random velocity component distributions (( pp = 1/3) and 
no figure rotation. 

In the extreme case where the reversion process is completed, n = N/2, 
the changes expressed by Eqs. (58)-(63) take the form: 

TV 

T rdm -> T rdm - ymf^)^] 2 ; (65) 

jV 

(T rdm )^ 0r -> (T rdm )^ r - ym[(^) w/2 f ; (66) 

(^rdm)« -> (^rdm)^ ~ ^MC^Wf 5 £ = S,t ] (67) 

TV 

^ s ^0 + -m[(!^y 2 ; (68) 

(^sysW -> + ym[(%) w / 2 f ; (69) 
TV 

(T sys )« -> + -m[(^),v /2 ] 2 ; £ = s,t ; (70) 

similarly, Eqs. (64) take the form: 

. (l/3)T rdm - {N/A)m[{v^) N/2 f 

Cu = (AT/0 , u x ^— ; £ = s,t ; (71a) 

C = (1/3)Trdm ■ (71b) 
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in any case, the anisotropy excess, C,u — (rr < 0, is counterbalanced by figure 
rotation. 

4.2 Axial velocity component reversion 

Performing the reversion process on a given fraction of particles, n/N, with 
respect to axial velocity components, implies the conversion of random (trans- 
lation) motion kinetic energy into systematic (translation) motion kinetic 
energy, as: 

1 2n 

A(T rdm ) rr = -A(T sys ) rr = --M(tv) 2 = -—nm[{v;) n ] 2 ; (72) 

where the remaining parameters are left unchanged. 

The reversion process implies the following energy changes: 

2n 

T rdm -> T rdm - — nm[(v;) n } 2 ; (73) 

2n 

(T rdm ) rr -> (T rdm ) rr - — nm[(lv) n ] 2 ; (74) 
2n 

T sys -> + — nm[(lv) n ] 2 ; (75) 
2n 

(T sys ) rr -> + — nm[(lv) n ] 2 ; (76) 

while the contributions from random motions along the x s and x t principal 
axis, (T rdm ) ss and (T rdm ) tt , remain unchanged. 

With the system being relaxed, £ = 1, in the case under discussion, the 
generalized and effective anisotropy parameters, ( pp and ( pp , coincide with 
their counterparts related to the usual tensor virial equations, and Eqs. (45a) 
and (46a) take the explicit form: 

f (l/3)T rdm 

U " T rdm - (2n/N)nmm n f ' U (&) 

, (l/3)T rdm - (2n/N)nm[(y;) n } 2 
W T rdm - {2n/N)nm[{v;) n f ' 

where the special case, n — 0, relates to the initial configuration, character- 
ized by isotropic random velocity component distributions (( pp = 1/3) and 
no figure rotation. 
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In the extreme case where the reversion process is completed, n = N/2, 
the changes expressed by Eqs. (73)-(76) take the form: 

N 
~2 

N 

(T rdm ) rr -> (T rdm ) rr - —m[(v;) N /2} 2 ; (79) 



Trdm T Tdm - — m[(v r ) N / 2 ] 2 ; (78) 



N 



r sys->0 + -ml(v;) N/2 y ; (80) 



; ; / 1 1 ( ■ , . i \ . ■ > 
" N ^ ,2 



(T sys )rr -»• + ym[(f r )w/ 2 ] | (81) 

similarly, Eqs. (77) take the form: 

(l/3)T rdm , . 

T T dm ~ {N/2)m[{v r ) N / 2 \ 

(l/3)T rdm - (jV/2)m[(^) jV/2 ] 2 

U T rdm - (iV/2)m[(w)^/ 2 ] 2 ' 1 J 

in any case, the anisotropy excess, Cu~Crr > 0, is counterbalanced by (imag- 
inary) figure rotation. 



4.3 Change of reference frame and imaginary rotation 

Performing the reversion process on a given fraction of particles, n/N, implies 
the conversion of random motion (translation) kinetic energy into systematic 
motion (translation) kinetic energy, with respect to the x r principal axis, 
according to Eqs. (72)- (76). The related kinetic-energy tensor component of 
the centre of inertia, by use of Eqs. (56) and (72) reads: 

1 2n 

(T c )rr = 2 M (^) 2 = j^nm[(v;) n } 2 = -A(T rdm ) rr ; (83) 

which, in the case under discussion, coincides with the kinetic energy of 
the centre of inertia, T c , being (T c ) ss = (Tc)tt — 0- I n t ne centre of inertia 
reference frame after the reversion process, the random motion kinetic-energy 
tensor component related to the (T r ' dm ) rr , by use of Eqs. (56) and (83), 

takes the expression: 

1 N n2 1 

( T rdm)rr = -771 ^ ~ "Cr = (T rdm ) rr - -Mi^f = (T rdm ) rr - {T C )rr 5 

1 i=l 1 

(84) 
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where the kinetic energy, Tq = (Tc) rr , is hidden by the change of reference 
frame (e.g., LL66, Chap. II, §8). 

Let the i-ih particle be at the distance, wf 1 = {[x^] 2 + [a^] 2 } 1 / 2 , from 
the x r principal axis, with velocity component, vf 1 . The imaginary angu- 
lar velocity (Caimmi 1996; C06; C07), iftW, can be defined in such a way 
the translational kinetic energy along the Xy clXIS IS counterbalanced by the 
imaginary rotational kinetic energy around the 

\m[v^] 2 + I m p«;«] 2 [iQ^f = ; (85) 

n?> = ; (86) 

where the index, i, is related to the i-th particle, and the factor, i, is the 
imaginary unit. In this view, the velocity components on the x r principal 
axis may be translated into imaginary tangential velocity components on the 
(0 x s x t ) principal plane, as: 

(hvJ 2 = (iw r n r ) 2 = -v 2 r ; (87) 

according to Eqs. (85) and (86). 

Let imaginary rotation around the Xy ClXIS be imparted to the particles 
where the reversion process has been occurred, and their counterparts with 
equal imaginary tangential velocity components, as prescribed by Eq. (87) 
particularized to the mean axial velocity component, v^, expressed by Eq. (56). 
The related increment in imaginary kinetic energy reads: 

1 2n 

A(T sys )^ r = -M(kv) 2 = —nm[(ivM 2 5 (88) 

and the combination of Eqs. (72), (83), (87), and (88) yields: 

^(T S y S )ci> r 4> r = — A(T sys ) rr . = — (T c ) rr = A(T rdm ) rr ; (89) 
the above results may be reduced to a single statement. 

Theorem 1. Given a R fluid with isotropic random velocity distribution 
and no figure rotation, let the axial velocity component reversion process 
be performed on a given fraction of particles, n/N, with respect to the x r 
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principal axis. Then turning to the centre of inertia reference frame with a 
kinetic energy loss, A(T sys ) rr = (Tc) rr , is equivalent to put the initial con- 
figuration into imaginary rotation around the x r principal axis, with square 
mean tangential velocity component, {iv^) 2 = — (v^) 2 , with a kinetic energy 
gain, A(T sys )^ r = -(T c ) rr . 

Accordingly, Eqs. (73)- (76) are replaced by the following: 

2n 

T sys -*0--™[(%) n ] 2 ; (90) 

(T sys W - - ^nm[(v^) n } 2 ; (91) 
1 2n 

(T sys ) u ^0--—nm[(y^) n } 2 ; (92) 

while the contributions from random motions remain unchanged, and the 
random velocity distribution remains isotropic (Cii = C22 = C33 — 1/3)- 

In the extreme case where the reversion process is complete, n = N/2, 
and the maximum amount of available imaginary rotation has been attained, 
the changes expressed by Eq. (90)-(92) take the form: 

^sys -> - ^-m[(v^) N/2 } 2 ; (93) 
TV 

-> - —m[(y^) N/2 } 2 ; (94) 

(T sys ) u -> - ^m[(v^) N/2 \ 2 ; (95) 

the above results may be reduced to a single statement. 

Theorem 2. Given a R fluid with isotropic random velocity distribution 
and no hgure rotation, let the axial velocity component reversion process be 
performed on a given fraction of particles, n/N , with respect to the x r princi- 
pal axis, and the reference frame changed into the centre of inertia reference 
frame. Then the resulting configuration with anisotropy excess, — ( rr > 0, 
Eqs. (77), is equivalent to the initial configuration with null anisotropy ex- 
cess, (u — ( rr = 0, and imaginary rotation around the x r principal axis, with 
square mean tangential velocity component, ijv^) 2 = —(v^) 2 . 
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4.4 Anisotropy excess and imaginary rotation 

Given a nonrotating (fl r = 0) isotropic = (22 — (33 = 1/3) configuration, 
let the tangential velocity component on the (0 x s x t ) principal plane be 
reversed in such a way Eqs. (65)-(71) hold. Let an equal amount of real and 
imaginary tangential velocity component on the (0 x s x t ) principal plane, be 
imparted to each particle for a total contribution equal to (2/3)/ r T rdm and 
— (2/3)/ r T rdm , respectively, where f r is a positive real number, which leaves 
the total energy unchanged. Let the reversion process be repeated for real 
tangential velocity components, to attain a null (real) figure rotation. The 
related changes, with respect to the initial configuration, are: 



rdm 



(l + -/ r ) T rdm ; (96) 

(T T dm)<t> r cl> r ^ (3 + 3^) ^ rdm ' (^) 
{T r dm)u — >■ - ^- + -fr^j ^rdm 5 (98) 

2 

-^sys ¥ g/r-^rdm i (99) 

2 

(^sysj^r^r ^ — -f r T rdm ; (100) 
1 2 

(T sys )i£ — »■ — --/ r T rdm ; (101) 
and the related anisotropy parameters read: 

[1 + (2/3)/ r ]T rdm 3 + 2/ r 1 ; 

C - (1/3)Trdm - 1 fi02b) 

U " [l + (2/3)/ r ]T rdm - 3 + 2/ r ' [i[)Zb) 

where the anisotropy excess, (u — (rr = /r/(3 + 2/ r ) > 0, is counterbalanced 
by imaginary rotation, according to an initial configuration with isotropic 
velocity component distribution and no figure rotation. As application of the 
above results, two significative examples shall be taken into consideration. 

First example. Nonrotating systems Rattened by anisotropic velocity com- 
ponent distribution (<7 n = cr 2 2 > ^33) ■ 
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Let tangential velocity components on the (OX1X2) principal plane be re- 
versed in a convenient fraction of particles, n/N, to yield a convenient amount 
of figure rotation together with isotropic velocity component distribution 
(a' u = a' 22 = 033), as sketched in Fig. 1. 

The combination of Eqs. (52), (53), and (55a) yields: 



while the mean square velocity components are left unchanged by the oc- 
currence of the reversion process. The substitution of Eqs. (103) and (104) 
into (37b) and (37c) shows the dependence of systematic and random motion 
tangential kinetic-energy tensor components on the reversion process. 

In the case under discussion (r = 3), the random velocity component 
distribution has to be isotropic after the reversion process, which makes 
Eqs. (60) and (63) reduce to: 



^nm[{v, r )S = l%Mimn? = f = y (4 - -33) I (107) 



which defines the configuration after the reversion process. 

Second example. Systems elongated by anisotropic velocity component 
distribution (<7 n = a 2 2 < 033) with no figure rotation. 

Let a convenient amount of real and imaginary figure rotation, and \v^ 3 , 
with respect to the X3 principal axis, be imparted to the system. Then the 
kinetic energy remains unchanged. Concerning real rotation, let the reversion 
process be performed on one half particles, leaving imaginary figure rotation 
together with isotropic velocity component distribution (a' n = <7 22 = 0-33), as 
sketched in Fig. 2. Accordingly, Eqs. (103) and (104) hold for imaginary tan- 
gential velocity components on the (Oil 12) principal plane and imaginary 
figure rotation around the x 3 principal axis. 




(103) 



(104) 




(105) 



(106) 
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Figure 1: After reversion of tangential velocity components in a convenient 
fraction of particles, with respect to the (Oa; 1 a: 2 ) principal plane, a configu- 
ration flattened by anisotropic velociiy2 component distribution (er n = a 22 > 
(T33) with no figure rotation (v^ = 0), left picture, is turned into a configu- 
ration flattened by figure rotation = (cr^ — o"^) 1 ^ 2 = (crfi — ^is) 1 ^ 2 ] with 
isotropic velocity component distribution (a' u = o' 22 = C33), right picture. 
The symbol, 13, denotes figure rotation around the related principal axis. 




Figure 2: After imparting a convenient amount of real and imaginary fig- 
ure rotation, and k^ 3 , with respect to the x 3 principal axis, to the 
system, and reversing real tangentiaP^elocity components on one half par- 
ticles, a configuration elongated by anisotropic velocity component distri- 
bution (o"n = 022 < 033) with no figure rotation (v^ = 0), left pic- 
ture, is turned into a configuration elongated by (imaginary) figure rotation 
[ii>0 3 = (cr^ — o^) 1 / 2 = (o"n — o"^) 1 ^ 2 ] with isotropic velocity component 
distribution (a' u = a' 22 = cr 33 ), right picture. The symbol, U, denotes figure 
rotation around the related principal axis. 



In the case under discussion (r = 3), an isotropic velocity component 
distribution after the reversion process implies the validity of Eqs. (105), 
(106), and (107), where the tangential velocity components are imaginary 
(au < 0-33), and the configuration after the reversion process is completely 
defined. 



4.5 Tangential velocity component reversion in the gen- 
eral case 

In the general case of anisotropic velocity component distribution (cr n 7^ 
^"22 7^ °"33) and figure rotation around each principal axis (fli ^ H 2 7^ 3 ^ 
0), with regard to the (Ox s x t ) principal plane, let the tangential velocity 
component reversion process be applied to one half particles in such a way 
no figure rotation around the x r principal axis occurs, Q r = 0. Keeping in 
mind Eqs. (50)- (53), the related energy changes read: 

(T rdm )^ r ^ r -> (T rdm )^ r + ^M(v^) 2 ; (108) 
(T sys W - (T sys )^ r - l -M{v^) 2 = ; (109) 



1 



N 



V<t, r = Mn/2 = T7$>2 5 ( 110 ) 
iV i=l 



A(T sys )^ r = -A(T rdm )^ r = -\m(v^) 2 = - l -MR 2 Gr {n r f ; (111) 
[{T r dm)<t>r4>r}u+[{Trdm) Wr w r ]u = ~Ma 2 u ; t — 1, 2 ; (112) 
(T rdm ) rr = ^Ma 2 r ; (113) 

where the rms velocity components, o 2 u and a 2 r , are related to the initial 
configuration. 

The changes in anisotropy paramenters, ( pp = a 2 p /a 2 , read: 

J- + A C« 

Crr^T^T- ; (H4b) 
1 + I\L,u 
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AC« = ^ ; £ = s,t ; (114c) 

where the rms velocity, a 2 , is related to the initial configuration. 

The application of the above procedure to the X\ and x 2 principal axes, 
makes the transition from an initial configuration with rms velocity compo- 
nents, o\i , a\ 2 i cr| 3 , and figure rotation around the principal axes, Qi, Q 2 , ^3, 
to a final configuration with rms velocity components, {c[i) 2 , {o'' 22 ) 2 , (°33) 2 > 
and figure rotation around the principal axes, 0, 0, f2 3 . The related energy 
changes read: 

(T rdm )n^(T rdm ) n + ^M(^) 2 ; (115) 
(T rdm ) 22 -> (T rdm )22 + ^M(v^) 2 ; (116) 
(T rdm ) 33 -> (T rdm ) 33 + \m [{v^) 2 + (v^) 2 ] ; (117) 

' ' ' r 'W7) 2 + (^) 2 l ; (us) 



(T rdm ) pp = \mo 2 w ■ p= 1,2,3 ; (119) 

and the related changes in anisotropy parameters are: 

c , C11 + (l/2)ACn . fl20 , 

Gl l + ACn + A C 22 ' (12 ° a) 

C22 + (1/2)AC 22 . ( 

C " l + ACn + A C 22 ' (12 ° b) 

C 33 + (l/2)[ACn + AC 22 ] . , , 

C33 ^+A Cll + A C 22 _ 2 (12 ° C) 

Adi = ^ ; AC 22 = ^ ; ; (I20d) 

the above results may be restricted to a single statement. 

Theorem 3. Given a J? .fluid, a convenient application of the tangential 
velocity component reversion process makes an adjoint configuration where 
figure rotation occurs around a single principal axis, that is a R3 fluid. 

Accordingly, the results valid for R3 fluids (C07) may be extended to the 
general case of R fluids. 
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5 Discussion 



As suggested in earlier attempts (Caimmi 1996; C06; C07), the equivalence 
between a variation in figure rotation and in anisotropy excess, may provide 
a useful tool for the description of collisionless fluids. The discussion here 
shall be focused on the spin parameter (Peebles 1969, 1971): 



where G is the gravitation constant, M the total mass, J the total angular 
momentum, and E the total energy. The above formulation includes four 
possibilities, namely (i) real rotation (J 2 > 0) and bound system (E < 0), 
which is the sole currently used in literature; (ii) imaginary rotation (J 2 < 0) 
and bound system (E < 0); (iii)real rotation (J 2 > 0) and unbound system 
(E > 0); (iv) imaginary rotation (J 2 < 0) and unbound system (E > 0). 
Accordingly, the spin parameter attains real values in cases (i) and (iv), and 
imaginary values in cases (ii) and (iii). 

In the light of the current model, oblate-like and prolate-like configura- 
tions belong to cases (i) and (ii) outlined above, while cases (iii) and (iv) rep- 
resent unbound structures for which the virial equations do not hold. Then 
the comparison with observations and/or computations, must be restricted 
to bound configurations. 

The spin parameter distribution is usually fitted using a lognormal dis- 
tribution (e.g., van den Bosh 1998; Gardner 2001; Ballin & Steinmetz 2005; 
Hernandez et al. 2007) or, in general, dependent on log A (e.g., Bett et al. 
2007), in dealing with real rotation. The inclusion of imaginary rotation 
would imply use of A 2 instead of A as independent variable, allowing for both 
positive (real rotation) and negative (imaginary rotation) values. 

The following procedure should be followed for calculating A 2 from obser- 
vations and/or computations: (1) determine the inertia tensor and the prin- 
cipal axes of inertia for an assigned matter distribution; (2) determine the 
potential-energy tensor; (3) determine the anisotropy parameters using the 
generalized virial equations; (4) perform the reversion process with respect 
to two principal axes of inertia to leave figure rotation around the third one 
(R3 fluid); (5) convert the anisotropy excess into (real or imaginary) figure 
rotation to obtain isotropic velocity component distribution (£ n = ( 2 2 — C33) j 
(6) evaluate the spin parameter; (7) act as already done for all the sample 



J 2 E 

G 2 M 5 
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objects; (8) determine the distribution of the spin parameter, P(A 2 ), with re- 
spect to the sample of adjoints configurations, where the velocity component 
distribution is isotropic. 

The related results could provide further insight on the formation and 
the evolution of large-scale collisionless fluids, such as galaxies and galaxy 
clusters. 

6 Conclusion 

A theory of collisionless fluids has been developed in a unified picture, where 
nonrotating (f^ = Q 2 = ^3 = 0) figures with isotropic (<7n = a 22 = C33) ran- 
dom velocity component distributions and rotating (f^ 7^ Q 2 7^ ^3) figures 
with anisotropic (er n ^ o 22 7^ 0-33) random velocity component distribu- 
tions, make adjoints configurations to the same system. R fluids have been 
defined as ideal, self-gravitating fluids satisfying the virial theorem assump- 
tions (e.g., LL66, Chap. II, § 10; C07), in presence of figure rotation around 
each principal axis of inertia. 

To this aim, mean and rms angular velocities and mean and rms tangential 
velocity components have been expressed, by weighting on the moment of 
inertia and the mass, respectively. The figure rotation has been defined as 
the mean angular velocity, weighted on the moment of inertia, with respect 
to a selected axis. 

The generalized tensor virial equations (Caimmi & Marmo 2005) have 
been formulated for R fluids and further attention has been devoted to ax- 
isymmetric configurations where, for selected coordinate axes, a variation in 
figure rotation has to be counterbalanced by a variation in anisotropy excess 
and vice versa. 

A microscopical analysis of systematic and random motions has been 
performed under a few general hypotheses, by reversing the sign of tangential 
or axial velocity components of an assigned fraction of particles, leaving the 
distribution function and other parameters unchanged (Meza 2002). 

The application of the reversion process to tangential velocity compo- 
nents, has been found to imply the conversion of random motion rotation 
kinetic energy into systematic motion rotation kinetic energy. The applica- 
tion of the reversion process to axial velocity components, has been found to 
imply the conversion of random motion translation kinetic energy into sys- 
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tematic motion translation kinetic energy, and the loss related to a change of 
reference frame has been expressed in terms of systematic motion (imaginary) 
rotation kinetic energy. 

A number of special situations have been investigated with further detail. 
It has been found that a R fluid always admits an adjoint configuration 
where figure rotation occurs around only one principal axis of inertia (R3 
fluid), which implies that all the results related to R3 fluids (Caimmi 2007) 
may be extended to R fluids. 

Finally, a procedure has been sketched for deriving the spin parameter 
distribution (including imaginary rotation) from a sample of observed or 
simulated large-scale collisionless fluids i.e. galaxies and galaxy clusters. 
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